Statistics of lowest excitations in two dimensional Gaussian spin glasses 
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A detailed investigation of lowest excitations in two-dimensional Gaussian spin glasses is pre- 
sented. We show the existence of a new zero-temperature exponent A describing the relative number 
of finite- volume excitations with respect to large-scale ones. This exponent yields the standard ther- 
mal exponent of droplet theory 9 through the relation, 6 = d(\ — 1). Our work provides a new way 
to measure the thermal exponent 6 without any assumption about the procedure to generate typical 
low-lying excitations. We find clear evidence that 8 < Qdw where Qdw is the thermal exponent 
obtained in domain-wall theory showing that MacMillan excitations are not typical. 
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Despite three decades of work in the field of spin glasses 
major issues such as their low-temperature behavior still 
remain unresolvedEJ. One of the main achievements has 
been a good understanding of mean-field theorycl which 
nevertheless does not include spatial effects which mani- 
fest as droplet excitations. Leaving aside the controversy 
whether replica symmetry breaking is or not a good de- 
scription of the spin-glass phaseu, it is reasonable to ex- 
pect that a phenomenological approach to the spatial 
structure of lowest excitations should satisfactorily ac- 
count fot, their low-temperature properties. McMillan 
proposecB that thermal properties in spin glasses are de- 
termined by the scaling behavior of the typical largest 
excitations present in the system. Therefore, he assumed 
that the energy cost of large scale excitations of length 
L scales like L e , 9 being the thermal exponent. Using 
domain-wall renormalization group ideas he also intro- 
duced a practical way to determine the leading energy 
cost of low-lying large-scale excitationsB. It consisted in 
measuring the energy defect of a domain-wall spanning 
the whole system obtained by computing the change of 
the ground state energy when switching from periodic to 
anti periodic boundary conditions in one direction. This 
idea has been further elaborated and extended to deal 
with equilibrium and dynamical properties of spin-glasses 
in a scenario nowadays referred as droplet model]. Sev- 
eral works have used McMillan's method Jo determine 
the value of 9 in two and three dimensionsQ. 

The purpose of this work is to show an alternative ap- 
proach to determine the low T behavior of spin glasses 
by studying the size and energy spectrum of the low- 
est excitations by introducing two exponents (A and 9') 
needed to fully characterize the zero-temperature fixed 
point. The entropic exponent A describes the probability 
to find a large-scale lowest excitation spanning the whole 
system, while the exponent 9' describes the system-size 
dependence of the energy cost of this excitation. 

The underlying theoretical background of our ap- 
proach is the following. To investigate the leading low- 
temperature behavior in spin glasses let us consider ex- 
pectation values for moments of the order parameter by 
keeping only the ground state and the first excitation. 



This approach was introduced irfl and can be extended to 
higher-order excitations to construct a low-temperature 
expansion for spin glasses.. This is described in detail 
in a separate publications. For simplicity, here we re- 
strict the analysis to the first linear expansion in T by 
keeping only the first excitation. If q = {a, r} denotes 
the overlap between two replicas (i.e. configurations of 
different systems with the same realization of quenched 
disorder), ..then the expectation value (q 2 ) can be written 
as followsQ, 
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where P(v, E) is the joint probability distribution to find 
a sample (among an ensemble of M s samples) where the 
lowest excitation has v spins overturned respect to the 
ground state (so the overlap between the ground and 
that excited state is q — 1 — 2v/V, V being the total 
volume of the system) and with energy cost or gap E. If 
v s and E(s) denote the volume and excitation energy of 
the lowest excitation for sample s then, 



(2) 



Using the Bayes theorem this joint probability distri- 
bution can be written as P(v,E) = g v P v (E) where 

J2v=i9v — Iq° dEP v (E) = 1, the last equality being 
valid for any v. g v is the probability to find a sample 
such that its lowest excitation has volume v and P V {E) 
is the conditioned probability for this excitation to have 
a gapL£qual to E. A simple low-temperature expansion 
of P)ErQ up to linear order in T yields, 
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showing that the leading behavior is determined by both 
the g v and the density of states at zero gap P„(0). In 
the droplet model it is generally assumed that typical 
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lowest excitations have average volume v = J2 V v 9v di- 
verging with V = L D and density weight at zero gap 
Py(0) ~ l/L 6 where 9 is the thermal exponent. In prin- 
ciple, a single exponent 9 describes the scaling behavior 
of large-scale excitations with volume v oc V which are 
supposed to determine the zero-temperature critical be- 
havior. Actually, assuming that gi/gv ~ 0(1), eq.(||) 
suggests that the contribution of small-scale excitations 
to the sum in the r.h.s would be negligible in the large 

V limit. Under these reasonable assumptions one ob- 
tains < q 2 > = 1 — cT/L e , where c is a non- universal 
stiffness constant related to the particular model. Our 
purpose here is to show that gi/gv diverges in the large 

V limit in such a way that the contribution of small-scale 
excitations to the r.h.s of (||) is as important as the con- 
tribution of the large-scale ones, thus leading to a new 
physical meaning of the thermal exponent 9. 
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FIG. 1. g(q) versus 1 - q for the PP (left panel) and FF 
case (right panel) for different lattice sizes L = 5 — 11 (PP) 
and L=6-16 (FF) from top to bottom. In both insets we plot 
the scaling function g(q)V x vs 1 — q with A = 0.7. 

Several numerical works have recently searched for 
low-lying .excitations in spin glasses using heuristic 
algorithmsEJ. But, to our knowledge, no study has ever 
presented exact results about the statistics of lowest exci- 
tations. In this paper we have exactly computed ground 
states and lowest excitations in two-dimensional Gaus- 
sian spin glasses defined by 

H = -} j Jij <7j a 3 (4) 

i<j 

where the cr^ are the spins (±1) and the Jij are quenched 
random variables extracted from a Gaussian distribution 
of zero mean and unit variance. These have been com- 
puted by using a transfer matrix method working in the 
spin basis. Representing each spins state by a weight and 
a graduation in the energy we can build explicitly the 
ground state by keeping the largest energy and next by 
iteration the first excitation and so on. The continuous 
values for the couplings assures that there is no accidental 
degeneracy in the system (apart from the trivial time- 
reversal symmetry a — > — a). Calculations have been 



done in systems with free boundary conditions in both 
directions (FF), periodic boundary conditions in both di- 
rections (PP) and free boundary conditions in one direc- 
tion but periodic in the other (FP). In all cases we find 
the same qualitative and quantitative results indicating 
that we are seeing the correct critical behavior. 

We found ground states and lowest excitations for sys- 
tems ranging from L = 4 up to L = 11 for PP and up 
to L = 16 for FP and FF. The number of samples is 
very large, typically 10 6 for all sizes. This requires a 
big amount of computational time and calculations were 
done in a PC cluster during several months. For each 
sample we have evaluated the volume of the excitation v 
(and hence the overlap q — 1 — 2v/V between the ground 
state and the first excitation) and the gap E. From these 
quantities we directly obtain g v and P V (E). In figure [j] 
we show g(q) — ^9v as function of q for different sizes 
in the PP and FF cases. We can clearly see that there 
are excitations of all possible sizes but the typical ones 
which dominate by far are single spin excitations. To 
have a rough idea of the number of rare samples giving 
large scale excitations let us say that nearly half of the 
total number of samples have one-spin lowest excitations, 
whereas less than 10% of the samples have lowest excita- 
tions with overlap q in the range — 0.5. This disparity 
increases systematically with size. A detailed analysis of 
the shape of g v reveals that it has a flat tail for large-scale 
excitations and a power law divergence for finite-volume 
excitations. The g v can be excellently fitted by the fol- 
lowing scaling form, 

- _ 1 (A B \ 
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This type of scaling, applied only to large-scale exci- 
tations, was proposed inu. Although we do not have a 
simple explanation for (|J) we emphasize how such ex- 
pression interpolates extremely well the whole spectrum 
of sizes matching the two volume sectors: a small-scale 
sector (finite v) g v ~ l/v x+1 and a large-scale sector 
(v/V finite) g v ~ 1/V +1 . This matching explains why 
there is a unique exponent A in (^|) which describes two 
completely different volume sectors. Note that although 
the g v is defined for discrete volumes, in the limit V 3> 1 
values of q are equally spaced by Aq — 2/V and the 
function g(q) = \g v becomes continuous if expressed as 
function of q instead of the integer variable v, 

^-2^( A+ (1-^TT) ■ (6) 

Although eq. (|5|) diverges for 5=1 apparently leading 
to a violation of the normalization condition for the g v , it 
must be emphasized that no excitation has q — 1 so there 
is a maximum cutoff value q* = 1 — 2/V corresponding 
to one-spin excitations. The normalization condition for 
the g(q) using (^) yields in the large V limit, 

fq *=i-2/v A- BIX B 
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implying A > 0. The divergent term (q — ► 1) in (g) shows 
that one-spin excitations dominate the whole spectrum. 
In fact, g(l) ~ 0(1) > g(V/2) ~ 1/V X+1 so the majority 
of excitations are finite- volume excitations. But the aver- 
age excitation volume v = Ylv=i V 9 V — AY 1 "* diverges 
in the V — > oo limit and differs from the typical excitation 
volume Vtyp ~ 1- This yields an independent measure- 
ment of the exponent A. By fitting the average volume 
to the expression v = C\ + C^V 1 ^^ or by measuring the 
ratio g(V/2)/g(l) ~ D 1 + D 2 V~ 1 ~ X we get an effective 
exponent A c ff = 0.7±0.05 as best fitting value. Although 
we find some systematic corrections in the value of A, our 
results are compatible with the general equation plus 
some systematic subleading corrections. We will see later 
how we can obtain a better estimate of A. 




FIG. 2. Gap distribution P{E) vs E for different lattice 
sizes in the PP case. In inset a) Scaling obtained from the 
ansatz (^) with 9' cff — —1.6. In inset b) we show the P V {E) 
for different excitation sizes (q = 0.5, q — 0) for a lattice size 
L — 10. Note that the distribution is independent of the size 
of the excitation. 

After having discussed the g v we show the results for 
the energy gap distribution P V (E). In figure || we show 
this distribution for the PP case. Similar results are ob- 
tained for the FF and FP cases. Quite remarkably, this 
distribution does not depend on the size v of the exci- 
tation, hence both large and small-scale excitations are 
described by the same gap distribution (inset b) of fig- 
ure ||). The normalized P V (E) has the following scaling 
behavior, 



P V {E) = L- 



(8) 



Since both the exponent 9' and the scaling function 
V are independent of v this implies that the gap prob- 
ability distribution P(E) = 'Yl lv>1 gvPv{E) satisfies the 
same scaling behavior (||) . In the inset a) of figure |^ we 
also show the best data collapse for P(E) obtained with 
an effective exponent 9' eS ~ — 1.6. A detailed study of 
the moments of P(E) computed for different values of L 
shows that there are also strong sub dominant corrections 
to the leading scaling (0) and that the value which gives 



the best data collapse turns out to be 9' = —1.7 ±0.1. We 
will argue below that, for Gaussian couplings, 9' < — 2 
and explain why finite-size effects are so big. 

Now we want to show how the exponents A and 9' 
combine to give the usual scaling exponent 9 describing 
the energy cost of typical thermal excitations in droplet 
theory. One of the most relevant results from the ansatz 
in ([)|) is that both small and large scale excitations con- 
tribute to low-temperature properties. In general, let us 
consider any expression (such as (^)) involving a sum over 
all possible volume excitations. Restricting the sum to 
the large-scale regime (v/V finite) the net contribution to 
such sum is proportional to VgvPv(fi) « L~ 6 ~ dx V(0). 
Coming back to (||) we note, using (||), that both small 
and large-scale excitations yield a contribution of the 
same order TL~ e ~ dx V(0). This is a consequence of the 
aforementioned fact that there is a single exponent A de- 
scribing the large and small sectors of the volume spec- 
trum (||). These considerations lead to the existence of an 
exponent 9 = 9' + dX determining the zero-tem perature 
critical behavior of the order parameter < q 2 >. Because 
of the systematic corrections in the values of the expo- 
nents A and 9' , in principle it is difficult to give an accu- 
rate estimate for the exponent 9. Nevertheless, despite 
of these strong corrections, the combination 9 = 9' + dX 
seems to be very stable (see figure (0J). Thus we com- 
puted, for each size, the combination 



v(L) 



(9) 



Since E(L) ~ L 6 ' and v{L) ~ L d ^~ x \ then the fit of 
A{L) gives a direct estimate of 9(L) = 6'(L) + dX(L). In 
Fig. we show the effective exponent form a fit of A(L) 
vs. L, with data in the range [L, ■ ■ ■ , L max — 11] for 
PPBC. We also show for comparison the effective expo- 
nent for 9dw obtained also with PPBC. Our best values 
for 9 is = —0.46(1) for the PPBC. This value is very close 
to the finite-temperature (Monte Carlo or transfer ma- 
trix) estimates 9 — —0.48(1)113 but certainly smaller than 
domain-wall calculations 9dw = —0.2850. Since our es- 
timate for 9 has been obtained with a new independent 
method it clearly shows that 9 ^ 9ow- 

Now we come back to the issue why in the large volume 
limit 9' should converge to a value < — d. The argumentO 
goes as follows. Consider the ground state and all pos- 
sible one-spin excitations. Because one-spin excitations 
are not necessarily the lowest excitations the statistics of 
the lowest one-spin excitations must yield a upper bound 
9[ for the value of 9' , i.e. 9' < 9[. The statistics of the 
lowest one-spin excitations is determined by the behav- 
ior of the ground-state local field distribution P(h) in the 
limit h — > 0. If P(h) is self-averaging and P(0) is finite 
in the large size limit then the statistics of the lowest 
excitations must be governed by the exponent 9' x = — d. 
Numerical resultsE3 in 2 dimensions show that the afore- 
mentioned conditions are clearly satisfied indicating that 
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9' < —2. Indeed, a determination of the finite-size cor- 
rections to the extreme statistics of lowest one-spin ex- 
citations shows that this asymptotic value can reached 
only for much larger sizes. The disagreement between 
our estimated effective value 6' cS ~ —1.7 in the scaling 
plot of figure 2 and the actually correct value 9' < — 2 
must be seen as a warning when extracting exponents 
from the average behavior of extreme values of statisti- 
cal distributions for small systems. 

If the power law distribution (^) properly interpolates 
the large and small volume sectors and then if we assume 
that 6' = —d holds in any finite dimension d, then the 
thermal exponent 9 is given by the interesting relation 
9 = d(\(d) — 1) indicating that X(di) = 1 corresponds 
to the lower critical dimension di. Above di the average 
lowest excitation volume is finite and A > 1. Our results 
may open the way to properly characterize the value of 
the thermal exponent 9 in three dimensions without need 
to use the MacMillan assumption about how to generate 
low-lying typical large scale excitations. Preliminary re- 
sults in three dimensionsEj reveal that (|J) describes well 
the data for small sizes. 
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FIG. 3. Domain-wall exponent (top) and 8 exponent (bot- 
tom) estimated as explained in the text and plotted as a func- 
tion of L. 

We have shown that a proper description of low- 
temperature properties in two-dimensional Gaussian spin 
glasses must be done in terms of two exponents: one 
(A) describing the relative number of finite-volume exci- 
tations compared to large scale ones, the other (9') de- 
scribing the typical energy cost of these lowest excitations 
whatever their size. Assuming that 9' = —d one can con- 
clude that the entropic exponent A fully characterizes the 
spin-glass phase. Although independent numerical esti- 
mates of 9' and A show strong finite-volume corrections, 
the thermal exponent 9 = 9' + dX is quite smooth with L 
giving 9 = —0.46 ± 0.01 unambiguously showing (given 
the high numerical precision for the value of 9dw) that 
9 < 9dw- To sum up, MacMillan excitations are not 
the typical low-lying excitations and our approach offers 



a new and independent way to estimate the thermal ex- 
ponent 9 without the need to generate typical low-lying 
excitations. Although we have focused our research in 
the two-dimensional Gaussian spin glass we believe that 
our conclusions remain valid for general dimensions be- 
yond d = 2. 
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